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Correlations in Many Electron Systems: 
Theory and Applications 

M. Tomaselli, T. Kuhl, D. Ursescu, and S. Fritzsche 



Abstract: In this contribution we present calculations performed for interacting electron 
systems within a non-perturbative formulation of the cluster theory. Extrapolation of the 
model to describe the time dependence of the interacting systems is feasible and planed. 
The theory is based on the unitary operator e lS (S is the correlation operator) formalism 
which, in this paper, is treated non perturbatively within many-particle correlations. The 
application of the derived equations to few-body systems is realized in terms of Generalized 
Linearization Approximations (GLA) and via the Cluster Factorization Theory (CFT). To 
check the reliability of the model we present two different applications. In the first we 
evaluate the transitions energies in Helium-, Lithium-, Beryllium-, and Boron-like Oxygen. 
The calculation aims to a precise determination of the satellite transitions which play an 
important role in plasma diagnostics. In a second we investigate a non-perturbative method 
to evaluate the charge radii of the Helium and Lithium isotopes by using the Isotopic Shift 
theory. We have found that our model leads naturally to components of e~ — e + pair in the 
two-electron wave functions of the Helium isotopes and three-electron wave functions of the 
Lithium isotopes. The possible connection of these terms to the QED leading diagrams is 
postulated. 

PACS Nos.: 31.10.+Z, 31.30.Gs, and 32.30.-r 
Resume : 

French version of abstract (supplied by CJP) 
[Traduit par la redaction] 

1. Introduction 

Deriving a non-perturbative and microscopic theory capable to describe the basic observable that char- 
acterize the dynamics of interacting electrons is a fundamental problem in the physics of atoms and 
ions. In general, one faces with two fundamental tasks, namely, the consideration of the correlation 
effects and the introduction of a cut-off parameter which, in order to obtain realistic and solvable sys- 
tems, reduces the dimensions of the model Equation of Motion (EoM). The introduction of correlation 
effects in many body systems via the e lS Unitary-Model Operator (UMO) goes back to the early work 
of Villars [1]. The idea is to introduce a wave operator S which maps zero-order reference wave func- 
tions (usually Hartree-Fock wave functions) to exact many body wave functions. Extended applications 
of the method in nuclear physics were shortly after performed by Shakin [2]. The e tS method came to 
quantum chemistry with the coupled cluster method proposed by Coester [3], and Kummel [4]. The 
coupled cluster Hamiltonian has been recently applied to the calculations of the electron affinities of al- 
kali atoms [5]. Studies of correlation effects in atomic systems based on the coupled cluster theory have 
been performed by Das et al. [6] . Recently [7,8] the e lS method was applied within nonperturbative ap- 
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proximations (Dynamic Correlation Model (DCM) and Boson Dynamic Correlation Model (BDCM)) 
to open shell nuclei. Applications of the method to open-shell electron systems were firstly applied to 
calculate the Hyperfine Splitting (HFS) constants of Lithium-like bismuth and uranium [10, 12]. The 
resulting non-perturbative and relativistic electron Dynamic Correlation Model (eDCM) was applied 
to calculate the effect produced by the electron and nucleon correlations into the isotopic shift theory 
IS. Calculations for lithium atoms were presented in [13]. Additionally the method finds application 
in the evaluation of dielectronic satellite-spectra of Lithium-like ions [14, 15, 16]. These are a useful 
tool for diagnostic of laser produced plasma. The ratio of various components of the satellite lines have 
been shown to be sensitive to density and temperature. 

We start by describing free electron systems with a relativistic shell model in which the wave func- 
tions are solution of the Dirac's equation. The model vacuum consists in paired electrons to fill major 
shells. The electrons in excess are considered as valence particles. The interaction between the elec- 
trons is responsible for exciting the valence electrons and for causing correlation effects in the closed 
shells. In additions to this polarization mechanism we have also the polarization of the continuum 
states. This polarization effects named Boiling of the Vacuum (BoV), have been already introduced 
in [10]. As in Ref. [7] we start by defining the basic operators of the model and by determining the 
relative EoM. The complex excitations modes are classified in terms of electron Configuration Mix- 
ing Wave Functions (eCMWFs). The eCMWFs form an orthogonal base of coupled clusters in which 
the Pauli principle between the different clusters is taken fully in consideration. Extrapolation of the 
non-perturbative cluster model to describe the time dependent electron-laser interaction is feasible and 
planed. 

In this contribution we present two applications of the non perturbative eDCM. The first involves 
the evaluation of the relativistic transition energies and wave functions for the Oxygen ions ranging 
from the Helium-like to the Boron-like. In the second application we study the dynamics of few- 
electron systems interacting with the excitation of the positron-continuum. The effect of this excitations 
is important in the determination of a non perturbative descriptions of the Mass Shift (MS) and Field- 
Shift (FS) which characterize the Isotopic Shift (IS) theory. 

2. Theory 

We start with a set of exact eigenstates {|^}} of the Dirac's Hamiltonian: 

hi = cdipl + ((3 - 1) + v nuc i(n) (1) 
which satisfies the dynamical equation 



In dealing with many electron systems one has to add the correlation effects caused by the two-body 
interactions: V(ij)c u\ an d V(ij)Brcit to the Hamiltonian of Eq. (1). Shell model calculation can be 
then performed to calculate transition energies between the different levels. Shell model calculations 
represent however an approximation in that one usually treats the effects of only few shells. The ne- 
glected shells serve to re-normalize the interaction in the shells considered. The re-normalization of the 
Hamiltonian is generally introduced via correlation operators. In UMO the effective Hamiltonian is cal- 
culated by introducing the correlations via the unitary e lS operator. By using only two body correlation 
we can derive: 



H\v)=E v \v) . 



(2) 



■ s *He iS * = 2<a|t|/3>4 + ^ (*a/j|«i2|* 7 «>4oJo 4 



a. 



'7 




(3) 
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where v\ 2 is the two body interaction and the ^ a p is the two particle correlated wave function: 

= e lS ^ af3 (4) 

However in dealing with complex atoms the (Si , i = 3 • • • n) correlations should also be consid- 
ered. The evaluation of these diagrams is, due to the exponentially increasing number of terms, difficult 
in a perturbation theory. 

We note that one way to overcome this problem is to work with e llySl+S2+s ^ operator on 

the Slater's determinant of the different states by keeping the n-body Hamiltonian uncorrected. 

After having performed the diagonalization of eigenvalue matrix obtained from the matrix elements 
of the n-body uncorrected Hamilton's operator, we can calculate the form of the effective Hamiltonian 
which, by now, includes correlation operators of complex order. 

The amplitudes of the correlated determinant are the calculated in the EoM method which is illus- 
trated in the following. 

If |0) denotes some physical vacuum and Oj, denotes the operator that creates the many-body 
eigenstate \u) such that Oj |0) = \u), O v \0) = 0, and H\0) = E o \0), then we have a set of EoM of the 
form 

ih lk\ 0) = [ ff ' ^°> = - Eo ^ = «>*Ol\0) . (5) 
In terms of the operators, the EoM can be written as 

[H,Ot\=u v Ot- (6) 
In Eq. (6) the Hamiltonian has the general second quantization form 

H = E a ^cic a + i J2 a ps-r( a P\v(r)\S 7)cj, cJcsg, 

= T + V int ' 

where T is the kinetic energy operator and Vi nt the interactions (Vcoui+^Breit). an d the d , c the general 
fermion operators. When they act on valence subspace, the and c creates and annihilates a valence 
electron, respectively. On the other hand, when they act on core subspace, the and c respectively 
annihilates and creates a hole state. Hence, the summation of the Greek subscripts leads to particle- 
particle, particle-hole, as well as hole-hole interactions. 

It is useful to determine the form of the central potential before the diagonalization of the model 
space is performed. This is because the matrix elements of the EoM can often be more easily calculated 
in a pre-diagonalization basis. 

^ {!&)}(— • • • , \r), . . .) is a complete set of basis vectors, then 

k)=5>}(%) = 5> 6l ,|&), (8) 



Ot = J2^0l. (9) 

6 

Using this last relation in Eq. (6), we obtain 

[h, oj] = ]T n br o\ (io) 

b 

where fl br = J2„ 

U r . Eq. (10) is the general form of EoM for the operator O v . 
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The coefficients Q(, r are simply the matrix elements of the Hamiltonian. To see this, we take the 
matrix element of both sides of Eq. (10) between the states (s\ and |0). Upon using the orthogonality 
between the basis vectors (i.e. (s|Oj|0) = S rs ), one obtains 



n sr = (a\H\r) - E S S 



(11) 



If the model space consists of a finite number, N, of basis vectors, then going from Eq. (10) back to 
Eq. (6) is equivalent to associate the systems of coupled equations given in Eq. (1 1) to the eigenvalues 
matrix equation given below: 



(O -El)x = 



(12) 



where O represents the (N x N) matrix O, 1 the ^-dimensional unit matrix, and x are the projections 
of the model space into the basic vectors. 

Equations (10) and (11) indicate that the complexity of solving Eq. (12) depends on the complexity 
of the model space, {\b)}, and the Hamiltonian, H. The following comparative review of the construc- 
tion of model spaces in different structure theories should give a glimpse on the scope of the problem. 

Let be the operator that creates n valence electrons outside the closed shells state $ ) : 



\m) = O^Xa^nh ■ ■ ■ in) = Yl a Jj $ °) = \ a "^hh ■ ■ ■ 3n) 



(13) 



In the simplest case where there is no closed shell excitation, the satisfies the EoM, Eq. (10) 

[H,Ol n (a m ;j 1 j 2 ---j n )] = ^2^mm'Ol n ,(a m r,j' 1 j 2 ---j' n ) , (14) 



with a m and a m > denoting the quantum numbers of the states |m) and \m'), respectively. 

The inert-core approximation would be good only if the valence-core interaction is very small. 
Hence, the applicability of the inert-core approximation is very limited as the interaction between 
valence and core electrons will generally excite the shell-model ground state of the core and create, 
in the process, the particle-hole (ph) pairs. Inclusion of the excitation mode due to lplh in the model 
space is known as the Tamm-Dancoff approximation (TDA) [17]. If one defines 



\m)TDA 



Al |0> 



TDA 



then Eq. (9) takes the form 



A in - 



3132 



Xjih a Ji b h®XZ 



(m')£)t 



(15) 



(16) 



The &] 2 creates a hole j 2 1 in |0)t_d.4 by destroying a core electron of j2 while a] i creates a valence 
electron of j\. The creates therefore a state of n + 1 particles and 1 hole (or p n+1 ft, 1 ). The x's 
are the configuration mixing coefficients and |0)t_d.4 denotes the physical vacuum of the TDA. In the 
literature one often chooses \0) T da = \HF), with \HF) being the Hartree-Fock ground state of the 
ion. In this latter case, 0^ m , = 1 in Eq. (16). 

It is also possible to use a physical vacuum that already contains ph pairs. In the literature, the 
method of random phase approximation (EPA) [17] has been introduced to study the full effects due to 
the pre-existence of lplh component in the physical vacuum. Hence, in RPA 



4> = E 



3132 



3l b 32 



+ X. 7, 7 , U 3 



(17) 
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and 

\m) R p A = RPA , (18) 

one can see that the term bj 2 aj 1 gives a null result if the physical vacuum | Orpa } does not contain pre- 
existing ph pairs. (In the literature, the coefficients Xhh an ^ X.hji are denoted by x™ j2 and —y 7 jlj 1 -) 
If the RPA is applied to closed-shell, then again Cft m , — 1 in Eq. (17). 

The introduction of the excitations of the vacuum in the above mentioned approximation is however 
complicated by the fact that the TDA and RPA vacua are different then the vacuum of the single 
particle operators. In addition simple calculations can be performed only by prediagonalizing the many 
body Hamiltonian in the TDA and RPA subspaces. The coupling to the additional valence particles 
can afterwards be accomplished by considering only few collective states and by neglecting the full 
treatment of the Pauling principle. In the following we show that these complications can be overcome 
by extending the EoM method to the field of non-linear equations. 

2.1. Polarization of the closed shells versus continuum vacuum excitations 

In the eDCM, the model space is expanded to include multiple ph excitations. This dynamic mech- 
anism includes either the excitations of closed electron shells or of positron-continuum states. More 
specifically [8], the eDCM states are classified according to the number of the valence electrons and 
of the electron particle-hole pair arising either from closed shells or from the positron-continuum. A 
state of N paired valence electrons and TV' particle-hole closed shells electrons ore" — e+ positron- 
continuum states is defined by 



\& N ^)=A\ NNI)J \0) (19) 
with 



^■(N,N')J 



X a N (J l .h...J N );J^ N {a N { J 1 J 2 ■ ■ ■ Jn); J) 

OIn(JiJ2---Jn) 

+ ^2 X a N + ll (J 1 .J 2 ...J N + 1 ,)-jA Ji N+v {a N+V (J 1 J2...JN+v)]J) 



a N + l'(JlJ2---J N + ll) 



(20) 



+ •••+ ^2 X aN+Nl (j 1 ...j N+N , ) . J A' < N+N ,(a N+N >(J 1 . . . J N +N'); J) 

a N + N<(Jl---JN+N>) 

where J denotes the total spin and the a's the other quantum numbers. The unprimed indices 1, . . . , n 
label the valence particle-particle pairs ( the valence bosons) and the primed indices 1', . . . , n' label 
the particle-hole pairs (the core electrons). The J/s denote the coupling of the pairs and the coupling 
of the different Ji is for simplicity omitted. The X's are projections of the model states to the basic 
vectors of Eq. (19). 

Within this definition the model space included either the excitation of the closed shells or the 
dynamics of continuum excitation which is taken into account through coupling the valence electron 
states to e~ — e+ states. The electron states defined in Eq. (19) are classified in terms of configuration 
mixing wave functions (eCMWFs) of increasing degrees of complexity (number of particle-hole or of 
e~ — e + pairs), see Ref. [7]. 

Since the different subspaces should be rotational invariant we introduce the coupling of the parti- 
cles and particle-holes in such a way that the first pair is coupled to angular momentum Jj , the second 
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to J 2, the two pairs are then coupled to J 3 and so on until all the pairs are coupled to the total angular 
momentum 7, e.g., 

A^ N (a N (J 1 J 2 ■ ■ ■ Jjv); J) = 



( [{ (44)^(44)^ } x h44V 3 ] Aa • • •) N_1 (<4-i4j j - 



(21) 



and 

^-N+l'( a N+l'(JlJ2 ■ ■ ■ Jn+1')J) = 

{( [{ (44) Jl (44V 2 } Al (44) j3 ] A2 • • • {4 n -An) Jn ) XN 
(«L+i'4„ +2 o J " +i '} J - 

Introduction of Eq. (20) into Eq. (10) gives the following equations of motion in the eDCM: 

[H, Aj v (aj V (JiJ 2 "-Ji V )J)]|0> 

= J2 Q p n p n A f N ((3 N (J 1 J 2 ---J N )J)\0) 

Pn 



(22) 



(23) 



(24) 



+ ^P N P N + 1 h 1 ^N+l'(PN+l'(JlJ2 ■ ■ ■ Jn+1')J)\0) 

0N + 1' 
+ ■■■ 

[H, A ] N+v {a N+v {J 1 J 2 ■ ■ ■ J N+V )J)]\Q) 

= ^■P N + I h 1 p N ^n(0n(JiJ2 ■ ■ ■ Jn)J)\0) 

Pn 

+ ^2 ^ l P N + 1 h 1 p N + 1 h 1 ^n+v{Pn+v{JiJi ■ ■ ■ Jn+1')J)\0) 

0N + 1' 

+ J]] ^p N + 1 h 1 p N + 2 h 2 A\f +v (0N+2>{JlJ2 ■ ■ ■ Jn+2>)J)\0) 

0N + 2' 
+ ■■■ 

where |0) is the shell-model state. Furthermore, we have used the notation p x h v for the indices of 
il to indicate the relevant xp — yh configuration. The additional commutator equations here are not 
given. In order to obtain eigenvalue equations we need to introduce a cut-off parameter: the GLA [7], 
which consists by applying the Wick's theorem to the A' n+2 ,(Pn + 2>( J\ J2 ■ ■ ■ Jn+2')J) terms and by 
neglecting the normal order. This linearization mechanism generates the additional terms that convert 
the commutator chain in the corresponding eigenvalue equation, as can be obtained by taking the 
expectation value of the linearized Eqs. (23) and 24) between the vacuum and the model states. 

Using the anticommutation relations and the Wick's algebra, one verifies easily that H can only 
connect states that differ by Iplh. The eigenvalue equation, Eq. (1 1), at the second-order linearization 
level is given by Eq. (25) where the subscripts referring to particle-hole configurations were not written 
explicitly but are understood. Note that in Eq. (25) ^l ap N pn p N+2 h 2 = Q a ,, p N+2 h2 p p N = 0. 



- E8 a f3 Q, a pi 

Q. a ip Cl a 'p> - 5 a >p> Sla'p" 

Sl a "P' Q, a "p" — E5 a »p» 



(25) 



The self-consistent method of solving Eq. (25) is given in detail in Ref. [8]. Here, we mention 
among others that in solving Eq. (25) the two-body interactions of H automatically generates nonlocal 
three-, four-interactions and so on. 
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The diagonalization of Eq. (25) can be performed only if one can calculate the many-body matrix 
elements. Calculations are feasible with the use of the Wick's algebra. However the number of terms 
to be evaluated increase exponentially and calculations are very slow. In this work, we perform calcu- 
lations by using the CFT of Ref. [7, 8, 11]. We believe that with the mastering of the essence of the 
CFT, matrix elements involving even more complex forms of operators can be easily deduced from the 
results obtained here. 

3. Transition energies in Oxygen ions 

The eDCM finds applications to the calculation of the transition energies of the Oxygen ions. In Table 1 
we give the energies for the Hydrogen-like Oxygen. The energies are calculate solving the Dirac's 
equation in a central Coulomb potential. For the Is i the calculated energy is compare wit the ionization 
energy of Ref. [18]. For the energies of the other levels no experimental energies are available. The 



Orbital 


Energy (eV) 


Ref. [18] 


Is 


-871.5080366004061 


871.41 


2s 


-217.9238060431288 




2p- 


-217.9234006166900 




2p 


-217.7378068079319 




3s 


-96.83498296247103 




3p- 


-96.83420091479068 




3p 


-96.77995294892040 




3d- 


-96.77911072701463 




3d 


-96.76111291825633 




4s 


-54.46202520546236 




4p- 


-54.46072932022532 




4p 


-54.43917924211600 




4d- 


-54.43737498444262 




4d 


-54.43087287667591 




4f- 


-54.42959174919736 




4f 


-54.42608486286169 




5s 


-34.85344986252919 





Table 1. Energies of the first 17 levels of the Hydrogen-like Oxygen. The minus sign designates the j — I — \ 
states. 

energies of the Helium-like Oxygen states are then obtained by solving Eq. (25). The indices (a, (3) 
are associated to a two electron states coupled to a good / quantum number. The energies of the first 
three J = + states, obtained by diagonalizing a matrix with 55 components, are given in Table 2. In 
the Table we give only three components of the 55 eCMWFs associated to the calculated spectroscopic 
factors. 



Spectroscopic factor 


Orbital 


Energy (eV) 


80% 


(lsilsi) u 


-1722.6 


87% 


(lsl2sl)° 


-1133.1 


79% 


(lsi3si)° 

2 2 


-1047.5 



Table 2. The first three levels of Helium-like O 6 with J=0 + and the associated spectroscopic factors. 

The energies of Lithium-like states are then obtained by solving Eq. (25). The indices (a, [3) are 
associated to a three electron states coupled to a good / quantum number. The energies of the first three 
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J = | states, obtained by diagonalizing a matrix with 350 components, are given in Table 3 together 
with the associated spectroscopic factors. 



Spectroscopic factor 


Orbital 




Energy (eV) 


92% 


{ls h {ls h 2p h Y] 


3 


-1949.2 


75 % 


(lsi (lsiSpi) 1, 

2 2 2 




-1849.1 


99% 


(lsi (lsi3p 3 ) lN 


3 


-1839.7 



Table 3. The first three levels of Lithium-like O 5 J=| and the associated spectroscopic factors. 

In order to calculate the transition energies of the Beryllium-like Oxygen we assume the first Is i 
shell full and we diagonalize Eq. (25) with the indices (a, (3) running over the unoccupied single 
particle states and the indices (a' , f3') over the lsi closed shell. The resulting energies for the three 
J = 1~ states obtained by diagonalizing a matrix of order 750, are given in Table 4 together with the 
relative spectroscopic factors. 



Spectroscopic factor 


Orbital 


Energy (eV) 


99% 


(2S12P1) 1 


-516.1 


98% 


((2si2si) l2 (2p 2 3 ls7 1 ) 1 ) 

2 2 2 ^ 


1 -434.6 


63% 


(2S12P1) 1 

2 2 


-434.3 



Table 4. The first three levels of Beryllium-like O J=l and the associated spectroscopic factors. 

In order to calculate the transition energies of the Boron-like Oxygen we assume the lsi shell 
full and we diagonalize Eq. (25) with the indices (a, [3) running over the unoccupied single particle 
states and the indices (ct',/3') over the closed shell. The resulting energies for the three J = + 
states obtained by diagonalizing a matrix of order 614, are given in Table 5 together with the relative 
spectroscopic factors. 



Spectroscopic factor 


Orbital 




Energy (eV) 


15 % 


(2si(lpi4si) 

2 2 2 




-5428.1 


37 % 


(2pi(3p|2pi) 1 (3pi 




-3897.3 


54% 


(2pi(2pi3pi) 1 (2pi 


2 S 7 1 )°)5 

2 


-3734.0 



Table 5. The first three levels of Boron-like O 3 J=§ + and the associated spectroscopic factors. 



4. Excitation of the positron continuum 
4.1. Non-linear realization of the IS theory 

The knowledge of the theoretical and experimental mass-dependence (MS) of selected atomic transi- 
tions and the theoretical calculations of the volume effects (FS) gives the possibility to have a deter- 
mination of the mean-square nuclear radii of short living isotopes [19]. Recent values for the nuclear 
charge radii of short-lived lithium and helium isotopes have been obtained from measurements per- 
formed at GSI, Vancouver [20], Argonne [21]. The measurements of the 2 2 S , 1 / 2 — > 3 2 5*i/2, of the 
2 2 S\/2 — > 2 2 P 3 / 2 , and 2 2 Si/ 2 — > 2 2 P!/ 2 transitions together with the recently performed calcula- 
tions [22] of the same transitions in lithium and helium atoms were in fact used to extract the difference 
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of the nuclear charge radii of the short-living isotopes from the charge radius of the stable isotope. In 
this paper we propose to reevaluate the MS and the FS in a non perturbative approximation based on 
the application of the eDCM. 

We start to calculate the energies of the lithium atoms by diagonalizing Eq. (25) in a base formed 
by three electrons in the (s,p,d) single particle states which interact with the BOV states formed by 
exciting the e~ — e+ continuum states. Results of this calculation for the 2s and 3s states are shown in 
table I. 



References 


method 


ls^2s energies in au 


Chung [23] 


FCPC 


-7.47805797(9) 


McKenzie and Drake [24] 


HR 


-7.478060326(10) 


Yan and Drake [22] 


HR 


-7.47806032310 (31) 


Tomaselli 


eDCM 


-7.478060733 


References 


method 


ls^3s energies in au 


Yan and Drake(02) [22] 


HR 


-7.354098422249(18) 


King [25] 


HR 


-7.354076 


Wang et al. [26] 




-7.5440980 


Tomaselli 


eDCM 


-7.35409801 



Table 6. Calculated energies of the ls 2 2s and the ls 2 3s in different models. 

According to Ref. [19] in order to evaluate the Mass Shift (MS) we have to add to the eigenvalue 
equation the additional term: Vj.Vj and to rescale energies and distance with the reduced mass of the 
electron. The matrix element of the Vj.Vj can be calculated as in Ref. [13] while the rescaling of the 
energies can be obtained by adding a i? n uci-^ clcctr ° n term to Eqs. (23) and (24) and to re-diagonalize 
the matrix given in Eq. (25). The correlations of the nucleus, which influence via the additional matrix 
elements given above, are in general approximated by a non relativistic perturbative calculations [22]. 

The FS term [28] factorize into a constant 

9tt7 

C = —m 3s - (Sr) 2s ] (26) 

where the term (5f) denotes the expectation value of the electron density at the nucleus multiplied 
by the isotopic variation of the charge radius. The polarizability of the nucleus, which influences the 
calculation of this constant, has been evaluated relative to the polarization of deuterium [27]. Since the 
FS is generally calculated in the point nucleus approximation, calculations performed within the DCM 
(nucleus) and the eDCM (electrons) correlation models could give a better insight in the FS calculation. 

Calculation of the IS for the isotopes of Lithium and Helium are under present calculation and will 
be reported soon. 

4.2. Transition energies in Lithium-like 235 U. 

The 2s-2p transition of Lithium-like 235 TJ is calculated in the eDCM. The result is given in Table 7 and 
compared with the QED calculation of Yerokhyn [29] and with the experimental result [30]. By using 
the resulting eCMWFS for the 2s i and 2pi we can calculate the hyperfine splitting (HF) of the two 
states. The calculation are performed by coupling the three electron wave functions to the ground state 
wave function of 235 U. For the nuclear ground state wave function we use a DCM which reproduce 
well within a large dimensional space the nuclear energies and moments of the 2/| valence neutron. 
Detailed calculation will be reported soon. 
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Transition 


2s-2p 


eDCM 


280.33 


Yerokhyn-Shabaev 


280.44(20) 


Experiment 


280.645(15) 
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Table 7. 2s-2p Transition in Lithium-like 35 U 



5. Conclusion and Outlook 

The transition energies of the Oxygen ions are calculated within a microscopic cluster model. The 
model is derived from the unitary operator model e lS which is used to generate the eCMWFs. The 
amplitudes of the model of the eCMWFs are calculated by using the EoM method. The modifications 
caused to the energy transitions by a time dependent laser beam can be simply evaluated by modifying 
the EoM. For few electron atoms like Helium and Lithium the EoM are extended to include the BOV 
excitations (excitation of the positron-continuum). The calculated eCMWFs together with the nuclear 
CMWFs of the different isotopes of Helium and Lithium allow a non-perturbative evaluation of the MS 
and FS of the IS theory. The influence of this new evaluation method on the charge radii of the Helium 
and Lithium isotopes is under present investigation. An open point in the presented calculation is the 
determination of the error of the calculated transition energies. In performing structure calculations we 
have used the single spinor energies obtained from the solution of the Dirac's equation (see Eq. (1). 
Better energies can be obtained by using the Harthee-Fock method. The approximation we have used 
gives to the calculated energies an error that can vary depending from the electron energy considered 
from 0. 1 to few percent. A better estimation of the errors could however be given, as suggested by 
Drake, by evaluating elementary excitation processes in light atoms like Hydrogen. For this purpose 
we are investigating the two photon transitions in Hydrogen. This would allow to establish a connection 
between the present non perturbative method and the QED perturbation theory. 
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